H = xp with interaction and the Riemann zeros 
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Starting from a quantized version of the classical Hamiltonian H = xp, we add a non local 
interaction which depends on two potentials. The model is solved exactly in terms of a Jost like 
function which is analytic in the complex upper half plane. This function vanishes, either on the real 
axis, corresponding to bound states, or below it, corresponding to resonances. We find potentials 
for which the resonances converge asymptotically toward the average position of the Riemann zeros. 
These potentials realize, at the quantum level, the semiclassical regularization of H = xp proposed 
by Berry and Keating. Furthermore, a linear superposition of them, obtained by the action of 
integer dilations, yields a Jost function whose real part vanishes at the Riemann zeros and whose 
imaginary part resembles the one of the zeta function. Our results suggest the existence of a 
quantum mechanical model where the Riemann zeros would make a point like spectrum embbeded 
in the continuum. The associated spectral interpretation would resolve the emission/absortion 
debate between Berry-Keating and Connes. Finally, we indicate how our results can be extended 
to the Dirichlet L-functions constructed with real characters. 

PACS numbers: 02.10.De, 05.45.Mt, 11. 10. Hi 



I. INTRODUCTION 

The Riemann hypothesis is considered the most im- 
portant problem in Analytic Number Theor y 1 ' 2 ' 3 ' 4 ' 5 . It 
states that the non trivial zeros of the classical zeta func- 
tion have real part equal to 1/2. Hilbert and Polya sug- 
gested long ago that the RH can be proved if one finds 
a self-adjoint linear operator whose eigenvalues are the 
Riemann zeros 6 -! 7 ^. The first indication of the adecuacy 
of this conjecture was probably the work by Selberg in 
the 1950s, who found a remarkable duality between the 
eigenvalues of the Laplacian acting on Riemann surfaces 
of constant negative curvature and the length spectrum 
of their geodesies^. Selberg trace formula, which estab- 
lishes that link, strongly resembles Riemann explicit for- 
mula. Another important hint came in 1973 from Mont- 
gomery's work who, assuming the RH, showed that the 
Riemann zeros are distributed according to the Gaussian 
Unitary Ensemble statistics of random matrix models 1 -^. 
Montgomery's results, were confirmed by the impressive 
numerical findings obtained by Odlyzko in the 1980's 11 . 
The next step in this direction was put forward by Berry 
who proposed the Quantum Chaos conjecture, according 
to which the Riemann zeros are the spectrum of a Hamil- 
tonian obtained by quantization of a classical chaotic 
Hamiltonian, whose periodic orbits are labelled by the 
prime numbers 1 ^. This suggestion was based on analo- 
gies between fluctuation formulae in Number Theory and 
Quantum Chaos 13 . Another interesting approaches to 
the RH are based on Statistical Mechanical idea o 14 ' 15 . 
The prime numbers has also been considered from a 
quantum mechanical viewpoint 1 - 6 -. 

Up to date, it is not known a Hamiltonian accom- 
plishing the Hilbert-Polya conjecture. Along these lines, 
Berry and Keating suggested in 1999 that the Id classical 
Hamiltonian H = xp is related to the Riemann zero s 17 ' 18 . 
This suggestion was based on a heuristic and semiclassi- 



cal analysis which yields, rather surprisingly, the average 
number of Riemann zeros up to a given height. Unfortu- 
nately, this encouraging result does not have a quantum 
counterpart. More explicitely, it is not known a quanti- 
zation of H = xp yielding the average, or exact, position 
of the Riemann zeros as eigenvalues. The Berry-Keating 
papers were inspired by an earlier one from Connes who 
tried to prove the RH in terms of the mathematical struc- 
tures known as adeles and p-adic numbers 1 ^. In order 
to illustrate the adelic approach, Connes introduced the 
Hamiltonian H = xp, using a different semiclassical reg- 
ularization. In Connes's approach the Riemann zeros 
appear as missing spectral lines in a continuum, which 
does not conform to the Berry-Keating's approach where 
the Riemann zeros appear as discrete spectra. Both ap- 
proaches are heuristic and semiclassical, therefore the ap- 
parent contradiction between them cannot be resolved 
until one derives a consistent quantum theory of H = xp, 
and its possible extensions. 

In reference 2 ^ we proposed a quantization of H = xp 
using an unexpected connection of this model to the one- 
body version of the so called Russian doll BCS model of 
superconductivit y 21 ' 22 ' 23 . The latter model was, in turn, 
motivated by previous papers on the Renormalization 
Group with limit cycle s 24 ' 25 ' 26 (see also 27 ' 28 ). The re- 
lation between H — xp and the Russian doll (RD) model 
is as follows. An eigenstate, with energy E, of a quantum 
version of the classical Hamiltonian H — xp, corresponds 
to a zero energy eigenstate of the RD Hamiltonian, where 
E becomes a coupling constant. Since the RD model 
is exactly solvable^ 3 -, so it is the H = xp model. The 
spectrum obtained in this way was shown to agree with 
Connes's picture of a continuum of eigenstates 20 . We 
also obtained the smooth part of the Riemann formula 
for the zeros, however this fact cannot be interpreted as 
missing states but rather as a blueshift of energy levels. 
A point like spectrum associated to the Riemann zeros 
was completely absent in this quantization of H = xp. 
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The final conclusion of 20 was the necessity to go beyond 
the H — xp model, in order to realize an spectral inter- 
pretation of the Riemann zeros. Some proposals were 
already made in that reference but the corresponding 
models could not be solved exactly. 

The cyclic Renormalization Group, and its realization 
in the field theory models of reference d 29 ' 30 ' 31 , is at the 
origin of LeClair's approach to the RH 32 . In this refer- 
ence the zeta function on the critical strip is related to 
the quantum statistical mechanics of non-relativistic, in- 
teracting fermionic gases in Id with a quasi-periodic two- 
body potential. This quasi-periodicity is reminiscent of 
the zero temperature cyclic RG of the quantum mechan- 
ical Hamiltonian o^ -, but the general framework of both 
works is different. The cyclic RG underlies several of the 
results of the present paper, but we shall not deal with 
it in the rest of the paper. 

The organization of the paper is as follows. In section 
II we review the Berry-Keating and Conncs scmiclassical 
approaches to H = xp. In section III we quantize this 
Hamiltonian, finding its self-adjoint extensions and their 
relation to the semiclassical approaches of section II. We 
also study the inverse Hamiltonian l/(xp) and its con- 
nection to the Russian doll model. In section IV we add 
an interaction to a quantized version of H = l/(xp), and 
solve the general model exactly, in terms of a Jost like 
function. Section V is devoted to the analiticity prop- 
erties of this Jost function. In section VI we study the 
potentials which exhibit some relation to the Riemann 
zeros. 



where A is the area of the allowed phase space region be- 
low the curve E = x p (see figs. lb, lc, Id). Table 1 collects 
the values oiAf(E) for the three types of regularizations. 



Type 


Regularization 


N{E) 


BK 


\x\ > l x , \p\ > l p 




C 


\x\ < A, \p\ < A 


flogA-£(log#-l) 


S 


l x < x < A 





Table 1.- Three different regularizations of H = xp and 
the corresponding number of semiclassical states in 
units h = 1. 

In the BK regularization, the number of semiclassical 
states agrees, quite remarkably, with the asymptotic limit 
of the smooth part of the formula giving the number of 
Riemann zeros whose imaginary part lies in the interval 

^^l( l0g |- 1 ) + ^ + --' E>>1 - {4) 

The exact formula for the number of zeros, Afjt(E), is due 
to Riemann, and contains also a fluctuation term which 
depends on the zeta function^, 



II. SEMICLASSICAL APPROACH 

The classical Berry-Keating-Conncs (BKC) 
Hamiltoniani2i±&±2 

H* = xp, (1) 

has classical trayectories given by the hyperbolas (see 
fig.la) 

x(t)=x e t , p(t)=p e- t . (2) 

The dynamics is unbounded, so one should not expect a 
discrete spectrum at the quantum level. In 1999 Berry 
and Keating on the one han d 17 ! 18 , and Connes on the 
other—, introduced two different types of regularizations 
of the model and made a semiclassical counting of states. 
Berry and Keating proposed the Planck cell in phase 
space: |x| > l x and \p\ > l p , with l x l p — 2irh, while 
Connes choosed |x| < A and \p\ < A, where A is a cut- 
off. In reference 2 ^ we considered a third regularization 
which combines the previous ones involving the position 
x, namely l x < x < A, making no assumption for the 
momenta p. The number, AT(E), of semiclassical states 
with an energy lying between and E is given by 



M R {E) = (N(E)) + Na(E) (5) 
(Af(E)) = ~ Imlog r Q + l -E^j - ?L log 7T + 1 

W fl (£) = = ilmlogcQ+iE) 

Based on this result, and analogies between formulae in 
Number Theory and Quantum Chaos, Berry and Keating 
suggested the existence of a classical chaotic hamiltonian 
whose quantization would give rise to the zeros as point 
like spectr a 17 ' 18 . They conjectured the properties of this 
classical Hamiltonian, which include the breaking of time 
reversal symmetry, which holds for |T]) , and the existence 
of primitive periodic orbits labelled by the prime num- 
bers. However, up to now, there is no a concrete proposal 
realizing all these conditions. 

On the other hand, Connes found that the number of 
semiclassical states diverges in the limit where the cutoff 
A goes to infinity, and that there is a finite size correction 
given by minus the average position of the Riemann zeros 
(see table 1). This result led to the missing spectral inter- 
pretation of the Riemann zeros, according to which there 
is a continuum of eigenstates (represented by the term 
— log A in Af(E)), where some of the states are missing, 
precisely the ones associated to the Riemann zeros. This 



3 




FIG. 1: la) a classical trayectory ([2}. The regions in shadow 
are the allowed phase space of the semiclassical regularizations 
of H — xp considered by: lb) Berry and Keating, lc) Connes 
and Id) Sierra. The values of the associated areas are given 
in table 1. 

interpretation, albeit appealing, has the trouble that the 
number of missing states changes linearly in E after scal- 
ing the cutoff A, and thus it is regularization dependent. 
As in the BK case, the C-regularization is not supported 
by a quantum mechanical version of H — xp, although it 
serves to illustrate, in a simple example, the main ideas 
underlying the Connes's adelic approach to the Riemann 
hypothesis. 

Finally, in the S-regularization the number of semiclas- 
sical states diverges as ^logA/Z^, suggesting a contin- 
uum spectrum, like in Connes's approach. But there is no 
a finite size correction to that formula, and consequently 
the possible connection to the Riemann zeros is lost. The 
main advantage of this regularization is that the Hamil- 
tonian |T]) can be consistently quantized yielding a spec- 
trum which coincides with the semiclassical result as we 
show below. 



III. QUANTIZATION OF xp AND l/(xp) 

A. The Hamiltonian Ho = xp 

In this section we construct a self-adjoint operator 
Hq, associated to Hq 1 = xp, which acts on the Hilbert 
space L 2 (a, b) of square integrable functions in the inter- 
val (a, b). Assuming that x > 0, there are four possi- 
ble intervals corresponding to the choices: a — 0, l x and 
b = A, oo, where l x and A were introduced above (we 
shall take l x = 1 and A = N > 1). Berry and Keat- 



ing defined the quantum Hamiltonian Hq as the normal 
ordered expression 



#o = -(xp + px), (6) 
where p = —ih d/dx. If x > 0, eq.© is equivalent to 

Hq = \fxp \fx = —ihy/x — \fx. (7) 

dx 

This is a symmetric operator acting in a certain domain 
of the Hilbert space L 2 (a, b), i££ 

(4>\H <I>) - (ffrW#) - in[ai>*{a)<t>{a) - bt/>*(b)(f>{b)] - 0, 

which is satisfied if both ^p(x) and <j>(x) vanish at the 
points a, b. By a theorem due to von Neumman, the 
symmetric operator Hq is also self-adjoint if its deficiency 
indices n± are equa^. These indices counts the number 
of solutions of the eq. 

Hl^ ± =±ih\ifj ± , (9) 

belonging to the domain of Hq (A > 0). If n = n + = 
?i_ > 0, there are infinitely many self-adjoint extensions 
of Hq parameterized by a unitary n x n matrix. The 
solutions of eq.© are 

i>±{x) = Cx~ 1/2TX , (10) 
whose norm in the interval (a, b) is, 



(^±)=±—(a^-b^). (11) 

The deficiency indices corresponding to the four intervals 
considered above are collected in table 2. 



Type 


(a, b) 


(n+,n_) 


Self-adjoint 


BK 


(l,oo) 


(1,0) 




C 


(0,iV) 


(0,1) 




S 




(1,1) 


V 


T 


(0,oo) 


(0,0) 


V 



Table 2.- Deficiency indices of Hq. The corresponding 

intervals are associated to the semiclassical 
regularizations of section II (i.e. BK, C, S). The last 
one, T, describes the case with no constraints on x 
except positivity (i.e. x > 0). 

The von Neumann theorem implies that the opera- 
tor Hq is essentially self-adjoint on the half line , = 
(0,co). This case was recently studied by Twamley and 
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Milburn, who defined a quantum Mellin transform using 
the eigenstates of H^. 

On the other hand, in the interval (1,N) the opera- 
tor Hq admits infinitely many self-adjoint extensions pa- 
rameterized by a phase e lS . This phase determines the 
boundary conditions of the functions belonging to the 
self- adjoint domain 

V(H , 9 ) = {ip,H iP G L 2 (l, N), e l V(l) = Vn^N)} ■ 

(12) 

The eigenfunctions of Hq, 

H ^ E ^E^, E , (13) 

are given byH 

M*) = -2i£m> ( 14 ) 

where C is a normalization constant. In the half line 
]K + there are no further restrictions on E, hence the 
spectrum of Hq is continuous and covers the whole real 
line ]R. In this case the normalization constant in (|14[l is 
choosen as C — \j \j2ixti which guarantees the standard 
normalization 



(0, oo) and (log a, log b), respectively. Under this map- 
ping, the wave function (fT4"|) corresponds to the plane 
wave e lqE , where x 1 ^ 2 is a measure factor. The spec- 
trum of Hq can therefore be understood in terms of the 
familiar spectrum of P. 

Returning to (flH|) . for the particular case where 6 = tt, 
one observes that the energy spectrum is symmetric 
around zero, i.e. if E n is an eigenenergy so is —E n . This 
result was obtained in reference 2 ^ working with the in- 
verse Hamiltonian 1/Hq. We shall next review that con- 
struction since it will be important in the sequel. 

B. The inverse Hamiltonian 1/Hq 

First, we start from the expression ([7]) and take the for- 
mal inverse, i.e. Hq 1 = x^ 1 ^ 2 p^ 1 x^ 1 ^ 2 . The operator 
is the one dimensional Green function with matrix 
elements (a;|p _1 |x') = ^ sign(x — x 1 ), where sign(x — x') 
is the sign function. The operator Hq 1 is defined in the 
interval (1, N) by the continuous matrix, 

H Q -\x,x') = ^ {gn{x - X '\ l<x,x'<N. (18) 
Its spectrum is found solving the Schrodingcr equation 



(V*')=C 2 J°°^x-« E - E 'V h = 6(E-E'). (15) dX 



, sign(a; — x 



'XX 



ip(x') = E- 1 ip(x), (19) 



In the case where Ho is defined in the interval (1, N), the 
boundary condition (112p yields the quantization condi- 
tion for E, namely 



N 



iE/h 



E n = 



2nh 
\ogN 



2n 



n € ]N- 



Hence the spectrum of Hq is discrete, with a level spac- 
ing decreasing for large values of N. The normalization 



constant of the wave function is now C = 
gives, 

(^E^eJ^C 2 I" *E. x -HE n -E n ,)/K 



which 



(17) 



The spectrum P^)l agrees with the semiclassical result 
given in table 1 for the S-regularization (recall that l x — 
l,A = N,H= 1). 

The existence of only two self-adjoint extensions of the 
operator Hq, in the positive real axis, should not be sur- 
prising, since they are intimately related to those of the 
momenta operator P = —ih-^, where q = \ogx. In- 
deed, the P operator defined on ]R, admits only two 
self-adjoint extensions in the g-intervals: (—00,00) and 
(a, b) (a, b finite), which correspond to the x-intervals: 



for the eigenvalue E , which must not be singular for 



.Hq to be invertible. Define a new wave function 



ip{x) 



(16) which satisfies 



iE 



JY 



2h 



Taking the derivative with respect to x yields 



x Tx m 



IT )"<•'■>• 



which is solved by 



= x l C lE /H =*" M*) = x l/2-iE/H > 



(20) 



— / dx sign(x — x') 4>{x) = x 4>{x). (21) 



(22) 



(23) 



with C = 1/VToglV as in (JTTJ). Eq. ^ fixes the func- 
tional form of i/j(x). To find the spectrum we impose (|2ip 
at one point, say x = 1, obtaining, 
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2irh 
logiV 



n6N (24) 



This spectrum coincides with (|16l) for 6* = 7r, so that the 
eigenenergies come in pairs {E n , —E n }, as corresponds to 
an hermitean antisymmetric operator. Including a BCS 
coupling in (fl~8|) . related to 9, yields the spectrum p6[) 20 . 

In summary, we have constructed in this section a 
quantum version of the classical Hamiltonian xp, as well 
as its inverse l/(xp), which agree with the semiclassical 
regularization l x < x < A. Both, the semiclassical reg- 
ularization, and the associated quantization, shows no 
trace of the Riemann zeros, which suggests that a pos- 
sible connection to them requires to go beyond the xp 
model. In the next section we shall take a further step 
in that direction. 



a(x) and b(x). If one of the potentials is constant, say 
b(x) — 1, then ([77)) becomes 



ATW) - - ^\ + < x l z^l (28) 



y/e{x)e(x') 



We shall denote M.i (resp. M.\) the model with Hamil- 
tonian (f27|) (reps. (|28|) ). These two models share many 
properties but they differ in some important instances. 
For example, (|27p is invariant under the transformation 



Vxe (i,N), 

(29) 

where the 2x2 matrix is an element of the «S7(2,R) 
group, 




IV. H = xp WITH INTERACTIONS 
A. Definition of the Hamiltonian 

The standard way to add an interaction to a free 
Hamiltonian Hq is to perturb it by a potential term, i.e. 



H = Hn 



V. 



(25) 



Instead of starting from the Hamiltonian Ho ^ we shall 
perturb the inverse Hamiltonian I/Hq 



a (3 
7 S 



£SZ(2,R)^a,/?, 7 ,«Se R, a£ - /3 7 = 1. 



while (1281) is invariant under the translations, 



a(x) — > a(a;) + a, 



ae R. 



(30) 



(31) 



Naturally, the spectra of Hi <x must be invariant under the 
corresponding symmetry transformations. The Hamilto- 
nian (|2"T|) can also be written as 



1 

H 



1 

Ho 



(26) 



so that H depends non linearly on V'. We have found 
more convenient to work with (|26|) , rather than with ([25)) 
but, of course, the two formulations must be related (we 
leave this issue for a later publication). 

The interacting Hamiltonian 1/H that we shall con- 
sider is given by, 



H^\x,x') = 



i sign(a; — x') + a(x)b(x') — b(x)a(x') 



2h 



v / e(x)e(x') 



... ( 2? ) 

where a(x) and b(x) are two real functions defined in the 
interval x € (1,N), and e(x) is a positive and monoton- 
ically increasing function. The BKC model corresponds 
to the choice e(x) — x, but it is equally easy to work with 
generic functions e(x), which also links the present model 
to the RD model, where e(x) gives the energy levels of 
electrons pairs. 

Hl^ x is an hermitean antisymmetric operator, and 
hence its spectrum is real and symmetric around zero. 
We shall assume, for the time being, that H^ 1 is in- 
vertible, a condition which depends on the potentials 



H 2 \x,x') = H Q 1 (x,x')+— [ij a (x)Mx') - Mx)Mx')} , 

(32) 



where 



a{x) b{x) 
ip a {x) = — ==, ipb{x) - 



^fe{x) 



y/eixj 



(33) 



This means that the interaction is given by a sort of 
proyection operator formed by the states ipa,b- In sec- 
tion VI we show that a particular choice of ty a .b provides 
a quantum version of the BK semiclassical regularization 
conditions. The Hamiltonians (f2"Tj) and (f2"B")) admit a dis- 
crete version analogue to the one considered ir>SP. The 
results we shall derive in the coming sections are also 
valid in this case. Furthermore, the connection with the 
RD model provides an interesting many-body generaliza- 
tion which will studied in a separate work. 



B. Solution of the Schrodinger equation 

The Schrodinger equation associated to (|2"?1) reads (in 
units of h = 1), 



6 



dx 



, sign(a; — x') + a(x)b(x') — b(x)a(x') 



'i y/e(x)e(x') 
which for the wave function 



<t>(x) 



^fe{x) 



= E-^ix) 
(34) 



(35) 



becomes 



iE rN 



e(x) 4>{x) 



(36) 



j l dx' (sign(x — x') + a(x)b(x') — b(x)a(x')) (j>(x'). 

This equation is the basis for the relation between the 
BKC and the RD models. Indeed, defining the RD 
Hamiltonian 



where 



dx a(x) 4>{x), B 



Eq. (|40j) is equivalent to 



iE 



TV 



dx b(x) 4>(x). 



(41) 



iEftx) + f x B-f x A=± (e(x) 4>{x)) (42) 



-if- J" dx 4>{x) + a x B - b x A = e(l) 0(1), (43) 
which are obtained from (|40|) by taking the derivative 



with respect to x, and setting x = 1 with a\ 
and bi — b(x ~ 1). Define the variable q 



dx 1 



e(x') ' 



a(x = 1) 



(44) 



such that 



Hrd 2 (x,x') = 



(37) 



e(x)S(x - a:') - % (sign(a; - x') + a{x)b{x') - b{x)a{x')) 



we see that (|36|) becomes the eigenequation of a zero 
energy eigenstate (j)(x), 



H 2 = E\4>) 



H 



RD 2 



(38) 



provided the coupling hp is related to the energy E by 



q G (0,L N ), 



L 



N 



N 



dx' 
e{x')' 



(45) 



In the BKC model, i.e. e(x) — x, one gets q — logx and 
Ln = logA^. For more general choices of e(x) we shall 
assume that Ln — > oo when N — > oo. In terms of the 
new function 



<p(x) = e{x) (j>(x), 



(46) 



hn = E. 



(39) 



Eqs.(|38|) and (|39j) establish a one-to-one correspondence 
between the energy spectrum of the Hamiltonian and 
the coupling constant spectrum of zero energy states of 
the Hamiltonian Hrjj 2 . In this regard, we shall mention 
the work by Khuri2&, based on a suggestion by Chadan 37 , 
where the Riemann zeros are related to the "coupling 
constant spectrum" of zero energy, s-wave, scattering 
problem for repulsive potentials in standard Quantum 
Mechanics. In that model the coupling constant, A, is 
related to the zeros, s n — 1/2 + ij n , by the quadratic 
equation A = s(s — 1). In our model however, the rela- 
tion between the coupling constant, hjm, and the energy 
E is linear (ea.(|39p). which is due to the fact that Ho de- 
pends linearly on d/dx, while in standard QM the kinetic 
term depends quadratically. 

We shall next solve eg. (|3"6"]) for generic values of a(x) 
and b(x). Later on, we shall impose additional con- 
straints on these functions so that the model makes sense 
in the limit N — * oo. First of all, let us write (|36|) as 



eq. p2|) turns into 



d S\ ~r, , da „ db , , 

d- q -* E )^ = T q B -d- q A > (47) 



where a, b and <j) are regarded as functions of q. If a 
and b are zero, the solution of (|47|) is the plane wave 
Ce iEq , with C a constant. Otherwise, C depends on q 
and satisfies 



dC 
dq 

whose solution is 



da db 
dq dq 



(48) 



C(q) = C+ P dq' e 
Jo 



- V / —B 



db 



dq' dq' 



-A 



(49) 



where C is an integration constant. This equation fixes 
the functional form of <p(q), namely 



r-N 



iE 

— j dx' sign(x— x') cj)(x')+a(x)B—b(x)A = e(x) <f>(x), 
2 Ji 

(40) 



<j){q) = Ce lEq + e lEq 



dq e 



i -iEq' 



da db ^ 



(50) 
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and in turn of by means of (|35|) and (|46l 



For e(x) = x, eq. (p)Tj) becomes 



dq' dq r 



(51) 



1 



„l/2-iE 



C 



(52) 

The term proportional to C in (|52p coincides with (|14p. 
The integration constants A, and C are related by 
eqs.([Iil and ([43jl. 



A = — 



dg a(g) </>(g), 



— y dq b(q) <f>(q), 



B 

* Jo 

•IT. <.LjV 

C = -— dq cj)(q) + ai B - hA. 

1 JO 

Plugging ([50]) into (|53|) yields 



(53) 



(1 + Ra ib ) A - R a<a B - R a C = 0, 
i? b , 6 A+(l-i? b , a ).B-i? b C = 0, 
(i2i,6 - 6i) A + (-i2i, a + oi) B - (iJi + 1) C = 
where 



(54) 



Rf(E) = 



iE 



iE 



dq /(g) e iE «, 



(55) 



We need in particular 



iEq 



t dg ; 



^3? 



Eq' 



The latter integral plays an important role in the sequel. 
Eqs. (|56[) can be proved using the change in the order of 
integration, 



dq I dq' = I dq' J dq. (59) 



o 



Plugging (O into (J3U) leads to 

(1 + S a , b ) A - S a , a B - R a (C + hA - aiB) = 0, 
S bib A+(l — S b , a ) B-R b (C + hA - aiB) = 0, (60) 
2e iEL N (ft b A _ + ( e iEL N + i)(C + hA - aiB) = 0. 

Observe that the term C + hA — a^B appears in all the 
equations. This parameter determines the asymptotic 
behaviour of the wave function at large values of q. In- 
deed, the value of C(q) at q = Lm (see eq. (|49| ) is given 
by 

Cv = C(q = L N ) = Coo + e~ lEL » (a N B-b N A), (61) 
where a at = a{x — N), 6jv = b(x — N), and 

Coo =C + A{h+2R b )-B{ ai +2R a ). (62) 
In the M.i model we shall impose 



lim a at = 0, lim 6jv 

N^oo N^oc 



0. 



(63) 



while for the A4i model, only the vanishing of a,oo will 
be required, since b(x) — 1. Under these conditions, Cjv 
converges asymptotically to Coo, as shown by (|6"Tj) . The 
latter parameter characterizes the asymptotic behaviour 
of the cigenfunctions, 



Ri(E) = l(e iELN ~l), (56) lim V(g)~4 = C 00 + ... (64) 




2 



Rij(E) = \h{l-e lEL »)-e lELN Rf{E), ^_ 

1 If Coo 7^ 0, the norm of tp(q) behaves as vXjv when 

Rf,g( E ) = Sf.g(E) — gi Rf(E), 

N — ► 00, and the wavefunction is delocalized, behaving 

, asymptotically like the "plane waves" of the non inter- 

wncrc 

acting model. On the contrary, if Coo = 0, the norm of 
ip(q) remains finite when N — > 00 for potentials decaying 
fi = f(x = 1), gx = g(x = 1), R f {E) = R f (-E), (57) sufficiently fast to infinity. An explicit expression for the 

norm of these localized states will be given at the end of 
and section V. The localization of the wave function -0(g) is 

due to an "interference" effect, which can dissapear by 
iE C Ln an infinitesimal change of the potentials. This suggest 

Sf ig (E) = — / dq f{q) g(q) (58) that the bound states of the model, if any, must be em- 

^ -'0 bedded in the continuum spectrum, unlike ordinary QM 

r n ^ f j ' 1 i\ -iEq' where the discrete and continuum spectra usually belong 

2 L 6 In ^ to separated regions. There are however exceptions to 
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this rule, as the class of von Neumann and Wigner oscil- 
lating potentials which have a unique bound state with 
positive energy embedded in the continuu m 33 ' 38 ' 39 ' 40 ' 41 . 

Continuing with our analysis, let us combine the third 
equation of ([50)) and ([6"2")l , obtaining 

C + b 1 A-a 1 B = -e lELN C OQ . (65) 

Plugging (|6"5"|) into (|6"0|) , yields the following linear system 
of eqs. 

Sw = 0, w* = (A, 5, Coo), (66) 

where 

/ 1 + S a , b Sa.a e* EL " R a \ 

S = S b , b 1 -S b>a e iBL »R b . (67) 
\ -2R b 2R a e tEL « + 1 J 

The existence of non trivial solutions of (|6"6"|) requires 

det S = T{E) + T{E) e lEhN = 0, (68) 

where 



3~{E) — 1 + S a , b — S b ^ a + S a ,a S byb — S a , b S b , a , (69) 
and 

T{E) — T{E) — 2R a (R b + R b S a , b — R a S b>b ) 

+2R b (R a — R a S b ^ a + R b S a , a ) ■ (70) 

To simplify J-(E), we use the equation 

SfAE) + S gJ (-E) = -2R f (E) Rg(—E), (71) 
which implies 



• If E is a solution of (|73j) . so is — E for generic po- 
tentials a and b, including complex functions. This 
is a consequence of the antisymmetry of H . 

• T{E) is invariant under the SL(2, ]R) transfor- 
mation ([29f . In fact, the terms S a , b — S b . a and 
S a ,a S b , b — S a , b S b , a are invariant separately. 

• We shall assume that the integrals defining Rf{E) 
and Sf y3 (E), converge for all values of E. At E = 
this implies 



Rf(0) = S f , g (0) = 0, /,S = a,6=^(0) = l, (74) 

Hence E = is not a solution of eq. ([73)1 and there- 
fore the Hamiltonian is non singular as as- 
sumed at the beginning of this section. 

Returning to the solution of eq. ([6"6")l . let us define the 
3-component vectors, 

Vl = (1 + S a , b , -S a ,a, e lELN R a ), 

v 2 = (S b , b , 1 - S b , a , e lEL » R b ), (75) 
v 3 = (-2R b , 2R a , e tEL » + 1), 

which are the rows of the matrix S, and are linearly de- 
pendent by (|68p . These vectors span a plane which can 
be characterized by its normal n. The vector w that that 
solves (|66p . must be proportional to n. If Vi and v 2 are 
non collinear, we can choose 

w = vi x v 2 , vi D" v 2 , (76) 

where 

vi x v 2 = (77) 
(-e iEL "{R b S a>a + R a (l-S bia )), 
e iEL »(R a S bib -R b (l + S afi )), T). 

We get in particular 



F{E) = T(-E). (72) 
The final form of the eigenenergies equation is 

T{E) + T{-E) e lELN = 0. (73) 

Before we analyze in detail the possible solutions of ([73)) 
we shall make some comments. 

• In the absence of interactions, i.e. a{x) = b{x) = 0, 
one gets T(E) = 1, VE, and then eq. ([73)) . for the 
case e(x) = x, reproduces eq. ([2"4")l . 



C 00 =T{E). (78) 

Hence the delocalized eigenstates, i.e. Coo ^ 0, satisfy 
that J-(E) 7^ 0, while the localized ones, i.e. Coo = 0, 
correspond to T{E) = 0. 

Otherwise, if vi and v 2 are collinear, the vector w can 
be choosen as 

w = vi x v 3 , if vi || v 2 , vi j| v 3 . (79) 

This case is exceptional since it requires the vanishing of 
the functions giving vi x v 2 . Another unlikely possibility 
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Bound states Scattering states 




* — » 



FIG. 2: Pictorial representation of the spectrum of the model. 
The bound states are the points where T(E) — 0, which are 
embedded in a continuum of scattering states. 



is that the three vectors Vj are collinear, in which case the 
vector w must belong to the plane orthogonal to vi^^. 

In summary, the most common situation we shall 
encounter is described by eqs. ([76|) . so that the local- 
ized/delocalized nature of the eigenfunctions is fully char- 
acterized by the vanishing/non vanishing of J~{E). The 
generic structure of the spectrum is depicted in figure [2] 
and table 3. 

In so far we have not used the reality of the potentials 
a and b, which guarantees the hermiticity of the Hamil- 
tonian (|27p . If they are real functions, then T(E) is a 
complex hermitean function, i.e. 

a{x),b{x)£^^F*{E)=F{-E*). (80) 
This equation follows from the identity 



C. Schrodinger equation for the Mi model 

The diagonalization of the Hamiltonian H\ (|28[) pro- 
ceeds along the same steps as for H 2 with the suitable 
changes. The main difference is that the phase factor 
e iEL N ar i ses m several expressions which one needs to ex- 
tract out and factorize conveniently. In particular, from 
([58)) one finds 

S a ,i = Ra, S lia = -e iEL » R a , S ltl = \{e lEL " - I), 
and then 



T = 1 + R a - -S a .a + e lEhN [Ra + ^S a .a + R a R a J , 

(83) 

so that eq. ([73|) becomes 

Fx (E) + e iELN fx (E) = 0, (84) 

where 

Tx{E) = l+2R a -S a , a , Fx(E) = l+2R a +S a . a +2R a R a . 

(85) 

Using (|7T|) for / = g = a, one finds the analogue of l|72[). 



S}jE) = S f », g *(-E*). 



(81) 



If E is real, eq.flHOJ) implies that T*{E) = T(-E) and 
hence the eigenvalue equation ([75]) . for delocalized states, 
can be written as 



HE) 
T{E)* 



for T(E) ± 



(82) 



The RHS of this equation describes the scattering phase 
shift produced by the interaction. For localized states, 
equation l|82p becomes singular, but eq. ([T3"]) is automat- 
ically satisfied since T(E) = J-{—E) = 0. All these re- 
sults means that 3-{E) plays, in our model, the role of 
a Jost function which determines completely the scatter- 
ing phases and bound states. However there are some 
important differences concerning their analytical proper- 
ties that we shall discuss below. 



Fx{E)=Fx{-E), 



so that 



Tx{E) 



AELn 



F X {-E) = Q. 



(86) 



(87) 



For the M.x model the constant is no longer related 
to the asymptotic value of Cjy (recall eqs. (|6"Tj) and l[S2"j)). 
Instead, it is replaced by the parameter Cx.oo which sat- 
isfies 



Cn — Cl.oo + e 



-iEL n 



a N B, 



Cx 

Cl,oo 



(88) 

eo = C^-e'^A, (89) 
iEL N (C + 2A-axB), 



where we have used eg. (1651) . The linear system (1601) for 
the constants A, B, Cx.oo turns into 



Eigenstate 


Coo 


HE) 


Eigencondition 


Delocalized 
Localized 


^0 
= 


= 


e F(E)* 

He) = o 



Table 3.- Classification of eigenstates of the M.2 model. 
For the Mx model, !F{E) is replaced by Tx{E). 



Siwi=0, w{ = (A, B, C liC 



where 



Si 



/ l + 2R a -S a , a e* EL »R a 
-1 1 - e *EL N 

2 2R n e iELN + 1 



(90) 



(91) 
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and whose determinant reproduces eq. (|87[) . 

det Si = Tx (E) + Tx (E) e lEhN . (92) 

Repeating the analysis made for the M.^ model we obtain 
the analogue of eq. (|78jl . 

Cx, 00 =T 1 {E), (93) 

so that Tx(E) controls the delocalized/localized charac- 
ter of the eigenfunctions. The unique exceptional case 
appears when Tx(E) = 0, R a (E) = -1 and e lELN = 1 
where, besides a localized solution with (A,B,Cx >00 ) = 
(1,1,0), there is the delocalized one (A,B,Cx t00 ) = 
(1, —1, —2). In summary, the generic eigenstates of the 
Mi model are described by table 3 with and T(E) 
being replaced by Cx )0O and T\{E). Most of the com- 
ments concerning eq. (|T3")) also apply to ([55)) . In addition 
we add: 

• The function Tx , unlike T , is not manifestly invari- 
ant under the transformation (|3ip . under which 



assume that these functions are well defined in the limit 
N — > oo, and call them T 00 (E) and Tx,oo{E). 

The reality of a(x) and b(x) guarantees the hermiticity 
of the Hamiltonians Hx,2, he. 

H\ 2 = Hx 2 a(x)* = a(x), b(x)* — b(x), Vx, 

_(97) 

and in turn that of the spectrum. Under these conditions, 
we shall prove the following important result: 



if ^ {E) = 0=> Im E < 0. (98) 

A similar statement holds for Tx, 00 (E). In other words, 
the zeros of the Jost functions, at N = oo, always lie ei- 
ther on the real axis, corresponding to localized states, or 
below it, corresponding to resonances. The proof of (|9"5|) 
is straightforward. Suppose that E satisfies eq. (|73l) . 
which implies that E is an eigenvalue of the Hamilto- 
nian Hi- Since the latter is hcrmitcan, E must be a real 
number, i.e. 



Tx ^ Tx + a 



JEL 



"(R a + I)-R a -I -a 2 (e lEL --l 



(94) 

However, eq. ([87|) is invariant under this change, as 
can be easily proved. In the large N limit we shall 
impose that oat — > 0, hence the symmetry (|3"Tj) will 
be fixed. 

• Eq.||7TD implies that 



if Tn(E) + Tn(—E) e iELN = 0=>ImE = 0. (99) 
Negating this implication yields, 



ME, L N if ImE ^0=>T N (E) + T N {-E) e ltjLN ^ 0. 

(100) 

Restricting to the case where ImE > and taking the 
limit TV — > oo (i.e. — > oo) in (|100|) one gets 



Tx(E) + Tx(-E) = 2(1 + R a (E))(l + R a (—E)). (95) 

If a(x) is real, the LHS of (|93|) gives the real part of 
Tx (E) , while the RHS is the norm squared of the 
function fx(E) (up to some constants) 

Re Tx (E) = \fx(E)\ 2 > 0, h{E) = l + R a {E). (96) 

Hence for the A^i model, the real part of Tx(E) is 
always positive. This property is analogue to the 
positivity of the imaginary part of the scattering 
amplitudes in some many body and QM models. 
On the contrary the real part of T can be positive 
or negative as we shall show in an example below. 

D. The zeros of the Jost functions 

The Jost functions T{E) and Tx(E) depend in general 
on the system size N. To make this dependence explicit, 
we shall denote them by Tn(E) and Tx : n{E). We shall 



ME, if ImE >0=>T aD (E)^ 0, (101) 

where the factor e lELN converges toward zero and thus 
cancells out the second term in (|100[) . Negating eq. (|101|) 
yields the desired statement (JHFJ- The case of Tx,oo{E) 
is similar. Repeating this argument in the case where 
ImE < does not give further information on the zeros 
of T(E). 

This property of the zeros of the Jost functions T(E) 
and Tx(E) is quite remarkable. In standard QM the ze- 
ros of the Jost function appear in the imaginary axis of 
the complex momentum plane (corresponding to bound 
states), or below the real axis (corresponding to reso- 
nances). Eq. (|9"8|) suggests a way to prove the Riemann 
hypothesis. Suppose, for a while, that T(E) were pro- 
portional to £(1/2 — iE). Hence, since the T(E) cannot 
have zeros with ImE > 0, the same property holds for 
£(1/2 — iE). The latter statement implies the RH. In 
section V, we shall see that C(l/2 — iE) does not have 
the correct analyticity properties to become a Jost func- 
tion of the model, but some modification of it may in 
principle. This issue will be considered in section VI. 
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FIG. 3: Plot of the real and imaginary parts of J-i(E), as 
given by eq. (|104[) . for the choices ai = 0.2,0.4,0.6,1. At 
ai = 1 the circle passes through the origin. 



E. Examples of Jost functions 

To illustrate the general solution of the Hamiltonians 
H\_2 we shall consider some simple models. 



Example 1: step potential in the Mi model 
Let us take 



a{x) — ai 8{x\ — x), 



1 < x < N, 



(102) 



where 9(x) is the Heaviside step function and 1 < x\ < 
N . R a and S a a are readily computed from eqs. ((55|) and 

(EE}, 



Ra = y (e lqiE - 1), S Q , a = y (e^ B - 1), (103) 

where qi = logxi. The associated Jost function J-\{E) 
follows from (l85l). 



Ti(E) = 1 



ai (2 - ax) 



(e lqiE - 1). 



(104) 



For each value of a\ ^ 0,2, the real and imaginary parts 
of J-\{E) describe a circle (see fig. [3]). At a\ = 1 the 
circle touches the origin, i.e. T\ = 0, at the energies 



EM = 

"2 



(2n 2 + 1)tt 



n 2 = 0,±l,. 



(107) 



In the thermodynamic limit, where Lm — > oo and qi 
is kept fixed, the spectrum of the Hamiltonian at a\ — 
1 consists of a continuum formed by the eigenenergies 
(|107p . and a discrete part formed by (| 1 05[1 (see fig. [2]). 
The explanation of these results is straightforward. The 
Hamiltonian H j~ , for the potential (|102p . has the block 
diagonal form 



H^(x,x') 



2y/e(x)e(x>) 







1 



\1 



«i 
«i 



1 — CLl 

1 — a\ 



a\ — 1 ai — 1 \ 
ai — 1 eti — 1 



-1 

1 

(108) 

where the vertical and horizontal lines separate the re- 
gions 1 < x < x\ and x\ < x < N, and 1 denotes a 
matrix with all entries equal to 1. At a\ = 1, the matrix 
(|108p splits into two commuting blocks whose structure 
is identical to that of Hq 1 in the corresponding inter- 
vals. Eqs. (|105|) and (|10T[) simply correspond to the non 
interacting eigenenergies in those regions. For ai ^ 1, the 
Hamiltonian H^ 1 is non diagonal and the eigenstates are 
plane waves in all the regions, with a phase factor e lqE 
and a discontinuity in the amplitude at x = x\. 

The function (|104[) is invariant under the transforma- 
tion 



ai <-> 2 — oi, 



(109) 



which maps a± = into ai = 2, and has ai = 1 as a fixed 
point. At ai = 2 one has T = 1, VE 1 , so that all the 
energy levels satisfy e lEhN = —1. One can check that 
the function is a plane wave e lqE , which changes its 
sign after crossing q = q±. More generally, one can define 
a unitary transformation which changes the sign of ^ for 
q > q±. Under this transformation the Hamiltonian for 
ai is mapped into that of 2 — ai, which explains eq. (|109p . 

Finally, let us look for the zeros of T\ (E) for generic 
values of oi, 



e n,E = _ x ^ E <J) = (2n i ±lK ; ni= o,±i 1 ... 

1 qi 

(105) 

describing an infinite number of bound states. The eigen- 
states corresponding to T\ ^ satisfy 



iEq x 



(ai - I) 2 + 1 
(ai 



1) 



1 



(110) 



For ai real, the RHS of (|110p has an absolute value 
greater than one, and therefore the imaginary part of 
E is a non positive number, 



iELf, 



(106) 



and their energy is 



F X {E) = and ai € R =J> Im£ < 0, 
in agreement with eq. 



(Ill) 
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Example 2: step potentials in the M2 model 
Let us take 



a(x) — a% 9(xi—x), b(x) = b\ 9{x2~x), 1 < x < N, 
where 1 < x\ < X2 < N. The S'-functions are given by 




Re F 



S a ,a = ^ (e^ E - 1), S M = ^-(e^ B -l), (113) 

„ ^ „ l FIG. 4: Plot of the real and imaginary parts of T(E) as 

S a , b = -y- (e iqiE - 1), S b , a = e 4(<?2 ~ 9l)B (e i<?l£; - 1) given by eq.CEIl for the choice c = 1/2 and qi/q 2 = 3/4 in 



where qi = logXi (i = 1,2). The Jost function is given 
by 



ea. (flTo|l . 



so that 



T{E) = l+c(c-l)(e lqiE -l)( e l ^- qi)E -1), 



ai&i 



2 

(H4) 

The conditions for this function to vanish, for real values 
of E, are 



?i(E) = 1 



This function vanishes at 



E + i[i 



(120) 



1 9i 



4ni ± 1 



C 2 ' 92 4n 2 + 1 ± 1 



which gives the eigenenergies 



, ni,n 2 € ]N, 



(115) 



.Fi(.E) = £ 



(1-f) 2 ' 



(121) 



in agreement with (|98|) , and it also has a pole at 



^m,n2 = 7^(4^1 ± 1) = T^(4n 2 + 1 ± 1). (116) 
2<?i 2q 2 



J-'i(E) = 00 ==> _E = — i/i, 



(122) 



which belongs to the lower half plane. The latter prop- 
FigHshows a particular example. To check that the zeros e rty is a general feature of the Jost functions, which con- 
of F{E) lie below the real axis for generic real values of s i st i n that their singularities always lie below the real 



<7i, q 2 and c one writes the equation T{E) = as 



axis. This result will be proved in the next section. 



i(qi-q\)E 



c(c- l)(e^ E - 1)' 



(117) 



If ImE > 0, the LHS of this equation is a complex num- 
ber with modulus less than one, in contradiction with the 
fact that the RHS has modulus greater than one. Hence 
one must have Im E < 0. 



Example 3: algebraic potential in the Mi model 
Let us choose 

a(x) = — , n > 0, 1 < x < 00. (118) 
x^ 

The R and S functions are given by 



Example 4'- algebraic potentials in the M2 model 
Let us choose 



a(x) = -^7-, b(x) = ^1,2 > 0, 1 < x < 00, 



with fii fi 2 - The £- functions are given by 



(123) 



S a .a 



a\ E 



4 E + ifii' 
a\b\iii E 



2(^i + /i 2 ) E + 
which yields 



Sb,b - 



E 



4 E + in2 



(124) 
E 



2(/xi + H2) E + i/j.2 ' 



Ra 



ai E 



2 E 



S a ,a 



4 E 



E —, (119) F(E) = l+(£-p) . f (125) 
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a\b\{ni - M2) 
2(mi + M2) 



(126) 



One can easily check that the zeros and the poles of this 
function always lie below the real axis. 

We have investigated other potentials to check ex- 
plicitely the property ([95)1 . For steps potentials, where 
the q intervals are related by rational fractions, eq. ([98|) 
follows from the Routh-Hurwitz theorem for the localiza- 
tion of the zeros of polynomials with real coefficients 4 ^. 
For more general potentials we have been able to check 
((98)1 numerically but not analytically. These results sug- 
gest the 7Wi j2 models may provide a huge class of com- 
plex functions with that interesting property. 



V. ANALITICITY PROPERTIES OF THE JOST 
FUNCTIONS 

In Quantum Mechanics the Jost function display 
analiticity properties which are a consequence of causal- 
ity. The close link between causality and analitic- 
ity is illustrated by the following theorem due to 
Titchmars h 43 ! 44 . 

Let /(g) be a generic complex function and f(E) its 
Fourier transform, 



f(E) 



dqf(q)e* E \ f(q)= I —f(E)e- iE ". 



(127) 

Assuming that f(E) is square integrable over the reals 
axis, 



dE\f(E)\ 2 <oo, 



(128) 



then, any of the following three statements implies the 
other two: 

• 1) /(g) = for q < 0. 

• 2) f(z) is analytic in the upper half plane, Im z > 0, 
and approaches f(x) almost everywhere as y — > 0. 
Further 



dx\f(x + iy)\ 2 <K, y > 0. (129) 



• 3) The real and imaginary parts of /, on the real 
axis, are the Hilbert transforms of each other, 



H[v], v = -H[u], f(x) = u(x) + iv(x), (130) 



H[g](x) = P [ * 

J-00 7T 



dy g{y) 



y-x 



(131) 



where P denotes the Cauchy principal value of the 
integral. 

Statement 1) is called causality, which in the present 
context means that the functions a(q) and 6(g) are zero 
for negative g-times, i.e. 



a(q) = b(q) = for q < 0. 



(132) 



In fact, the variable q is always non negative by ea.(|45p. 
so that (|132p must be understood as an extension of the 
definition of a(q) and b(q) for negative values of q. The 
main consequence of causality are the dispersion relations 
(|130p . which play a central role in the scattering theory 
in Quantum Mechanics, and other fields of Physics. 

For the M2 model to be well defined in the limit N — » 
00, we shall impose that / = a,b are square integrable 
functions, i.e. 



dq\f{qf 



dE 
2^ 



\f(E)\ 2 < 



(133) 



Hence by the Titchmarsh theorem, a(E) and b(E) are 
analytic functions in the upper half plane and satisfy 
eq. (|130p . which can be combined into 



f(E) 



dt f(t) 



in t-E 



f = a,b. 



(134) 



These properties, in turn, imply that Sf i9 and T are 
analytic functions. 



A. Analyticity of Sf lB 

Consider the .S-function defined in 

N -> 00, 



in the limit 



iE f°° 

SfAE) = Y J dqf(q)g(q) (135) 

-=T / dqf(q)e iE « / dq 1 g(j)e^ E \ 
1 Jo Jo 

where / and g are causal functions, in the sense of (|132|) . 
and square normalizable. Replacing / and g by their 
Fourier transforms one arrives at 



s f , 9 (E) = f r ' ^f(t)g(-t) 



(136) 



— 00 
2 roo 
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Using (|134p for g yields, 



SfAE) = -~f{E) g(-E) + — — t _ E , 



or equivalently 



4 S Lg (E) 
E 

-- Ef(E) g(-E) 



(137) 

(138) 
dt t f(t) g(-t) 

7T t-E 



which shows that Sf } g(E) is the sum of the function 
k(E) = tf{E) g(—E) and its Hilbert transform, i.e. 



Sf,g = k — iH[k]. 



(139) 



If k(t) belongs to the space L p = L p (—oo,oo), then Tt[k] 
is defined and belongs to L p for p > 1—. In this case the 
Hilbert transform of s/ i9 satisfies eq. (ll34[) . Moreover, if 
SfAE) is square normalizable then, by the Titchmarsh 
theorem, it will be an analytic function in the upper half 
plane, i.e. 



If tf{t)g{-t) e L 2 =*> Sf, g : analytic in (D+ (140) 

Apparently the normalizability of / and g does not guar- 
antee that of s/, 9 , but it all the examples we have ana- 
lyzed that is the case. 



B. Analyticity of Ti(E) 
A consequence of eg. (11401) is 



if t\a(t)\ 2 e L 2 



ME) - 1 

E 



analytic in (TJ + . (141) 



To prove ([Mil) , write T\{E) as (recall eqs. ([85]) and (J55 



ME) = l + iE a(E) - S a a (E). 



(142) 



Recall that a is analytic in (fj , as well as S a At)/t, pro- 
vided that t \a(t)\ 2 G L p> i. Hence (|141|) follows. We 
also expect that under appropiate conditions on the po- 
tentials, the combination {T{E) — 1)/E will be analytic 
in the upper half plane. 

We proved in section IV that the real part of T\ (E) is 
always positive and equal to the square of the function 
h{E) (see m), 



h(E) = l + ^-a(E). 



(143) 



Using eq. (|142p and the analyticity of a one can prove 
that the imaginary part of J-\ (E) is given by the Hilbert 
transform of \fi(E)\ 2 , i.e. 



ME) = \h(E)\ 2 -iP [ - 



dt \h(t)f 



t-E 



(144) 



However, the converse is not true. The reason is that 
T\ (E) does not in general converge toward zero as Im E 
goes to infinity. This can be simply illustrated by the 
non-interacting case where J-\(E) = 1. 

Finally, we shall give the expression of the localized 
eigenstates of the M.2 model in the limit Ljy — > oo, i.e. 



(^e\iPe) = 




(145) 




where 



^1,9 = - 2 ( Sg,qf + Sf. qg + — (SgJ - Sf,g)\ ~J rf? /ff. 

(146) 

for /, g — a, b. Using the Fourier transforms of these 
functions one can write (|146[) as 




(148) 



where d /dt only acts on the function f(t)*. The re- 
sults obtained on this section can be given a more formal 
treatment using the Theory of Hardy spaces^, but we 
leave this more mathematical matters for another work. 
In connection to the previous discussions we would like 
to mention the work by Burnol, who has emphasized the 
importance that causality in scattering theory may play 
in the proof of the Rlemann hpothesi a 47 ' 48 . 
Let us consider again some examples of potentials in- 
spired by the previous results. 



Example 5: Other algebraic potentials in the A4i model 
Let us choose fi(E) as 

2M 

h{E)=C 1 +C 2 \{ a " + d + C 2 = l, (149) 



a n — iE 
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where a n > (n = 1, . . . ,2M) to avoid the poles of 
fi(E) in the upper half plane. The condition C\ +C 2 = 1 
guarantees that fx (0) = 1 . From (|143p one has 



a(E) 



'2M 

n 



iE 



2C 2 



\n— 1 



(150) 



The quantity in parenthesis is an analytic function in (C + 
which behaves as 1/E for |E| >> 1. The associated po- 
tential a(q) can be computed from the Fourier transform 
ofa(E) 



2AI 



2M 



a( 9 )=4C 2 ^e-^ J] 



(151) 



n=l 



and consists in the superposition of decaying exponen- 
tials. In terms of the x variable (s(x) = x) the decay is 
algebraic. Using ()149[) one can show that 



Fx(E) = 1-2C 1 + 2C 1 / x (£), 
and in particular 



(152) 




Re F 



FIG. 5: Plot of the real and imaginary parts of T(E) as- 
sociated to the potentials (|157l) and (|159|) . E varies in the 
interval (1.91,2.05). The curve passes through the origin at 
E = 1.95634. 



where a n ,(3 m > 0. The S'-functions can be easily com- 
puted using 



S 



h,h 



Ml 



E 



2(/ii + E + ifii 



fi = e->»« (i = 1,2). 

(158) 

Fig. [5] displays an example with N a = 2 and Nj, — 1 and 
the following values of the parameters in (I157j) 



C x = C 2 = - Tx{E) = h(E). (153) 



The zeros of T\ (E) are given by 



2M 

n 



a„ + iE _ 1 - 2d + 2d 2 
a„ - iE ~ 2Ci(l-Ci) 



(154) 



For Ci = 1/2 there are 2M real solutions which appear 
in pairs {-E, — -E}, while for other values the solutions are 
complex and satisfy Im E < in agreement with ([55)) . 
When M = 1 and C\ = 1/2, the potential (|151[) becomes 

a(g) = 2 Q2 + Q!l {e- aiq ~ e~ a2q ), (155) 
a2 — ai 

and the solutions of (| 1 54[) are the pair of energies 



E = ±^axa2- 



(156) 



Example 6: Other algebraic potentials in the M2 model 

A general choice of the potentials a and b for the 
model is given by 



a n e 



6(g) = £ 6„ 



(157) 



ai = 1, a 2 = 4, ft = 0.0360157, (159) 
01 = 7.22928, a 2 = -7.03245, 61 = 1. 

At E = ±1.95634 the Jost function vanishes. This exam- 
ple is a perturbation of the potential (|15ip with ax = 1 
and c«2 = 4, which has bound states at E = ±2. The 
most interesting feature of this example is that Re F(E) 
becomes negative in a small neighbour of the origin. This 
has been possible by the addition of the b potential. No- 
tice that ax , 02 do not balance exactly as in eq. (|155|) . 



VI. THE RIEMANN ZETA FUNCTION AND 
THE JOST FUNCTION 

There are two main physical approaches to the Rie- 
mann zeros, either as a bound state problem, or as a 
scattering problem. In the former approach one looks 
for a Hamiltonian whose point like spectrum is given by 
the Riemann zeros, while in the latter the Ricmann zeta 
function gives the scattering amplitude of a physical sys- 
tem, whose properties reflect in some way or another the 
existence of the zeros. Both approaches would naturally 
converge if the Riemann zeros were the zeros of a Jost 
function as suggested above. 

The scattering approach was pionered by Faddeev 
and Pavlov in 1975, and has been followed by many 
author a 49 i 50 ' 51 i 52 . An important result is that the phase 
of C(l + it) is related to the scattering phase shift of 
a particle moving on a surface with constant negative 
curvature. The chaotic nature of that phase is a well 
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known feature. Along this line of thoughts, Bhaduri, 
Khare and Law (BKL) maded in 1994 an analogy be- 
tween resonant quantum scattering amplitudes and the 
Argand diagram of the zeta function C(l/2 — it), where 
the real part of £ (along the x-axis) is plotted against 
the imaginary part (y-axis)^. The diagram consists of 
an infinite series of closed loops passing through the ori- 
gin every time £(1/2 — it) vanishes (see fig. [B]). This loop 
structure is similar to the Argand plots of partial wave 
amplitudes of some physical models with the two axis 
being interchanged. However the analogy is flawed since 
the real part of ((1/2 — it) is negative in small regions 
of t, a circumstance which never occurs in those physical 
systems. 

In fact, the loop structure of the models proposed by 
BKL is identical, up to a scale factor of 2, to the model 
of example 1 (see fig. ([3])), where the loops represent- 
ing !Fi(E), for a\ = 1, are circles of radius 1/2, centered 
at x = 1/2. For general models of type Mi, the loops 
are not circles but the real part of T\ (E) is always posi- 
tive (see eq.(jnU)), and therefore they can never represent 
£(1/2 — iE). Incidentally, this constraint does not apply 
to the models of type M2, where Re !F(E) may become 
negative, as in the example 6 (see fig. [5]). This suggests 
that ((1/2 — iE) could be perhaps the Jost function T(E) 
of a M2 model for a particular choice of a and b. How- 
ever, ((1/2 — iE) has a pole in the upper half plane at 
E = i/2, while J-(E) is always analytic in that region. 
The solution of this problem consists in moving the pole 
to the lower half plane defining the function 




Re z 




FIG. 6: Left: real and imaginary parts of ((s) with s — 1/2- 
iE and E £ (0,50). Right: same as before for (h(s) = (s - 

i)C(«)/«. 



o 2i0(t) 



2 / 



it "i 
2 / 



(163) 



which is taken to be continuous across the Riemann zeros. 
This angle gives the smooth part of the Riemann formula 



i.e. 



7T 



(164) 



Ch(s) 



1 



C(s) 



(160) 



which has a pole at E = —i/2 for s = 1/2 — iE. This 
function was already considered by Hardy, and it is dis- 
cussed in detail by Burnol in his approach to the R H 47 i 48 . 
FigJB] shows the Argand plot of (#-(1/2 — it). In the rest 
of this section we shall explore the possibility that the 
zeta function, or some other related function, could be 
realized as a Jost function. 



As noticed by BKL, the loop structure depicted in fig. 
[6] shows that the zeros of Z(t) are near to the points 
where ((1/2 + it) is purely imaginary, i.e. 9(t) = tt(ji + 
1/2)^. This observation suggests an approximation to 
the Riemann zeros 



cos 9(t) = 



2 I 



= 0, 



(165) 



A. The Bessel potentials and the smooth part of 
the Riemann formula 



The zeta function satisfies the well know functional 



equation, 



1 



7T 



C( it) - 

l 2 ' r(i 



/ 1 i 2 ^ (( 



f) -2 



it) 



For t real, one defines 



(,{1/2 + it) = Z(t) e 



i9(t) 



(161) 



(162) 



where Z(t) is the Riemann-Siegel zeta function, which is 
even and real, and 9(t) is a phase angle given by 



which works within a 3 % of error (see fig. [7]) and it 
is essentially the same as the smooth approximation dis- 
cussed in section II. Condition (|165|) was also obtained by 
Berry from the first term in his approximate formula 12 . 
BKL related 9(t) to the scattering phase shift of a non 
relativistic particle moving in an inverted harmonic os- 
cillator, H = p 2 — x 2 , a problem which is related to the 
H — xp by a canonical transformation. 

In this section we shall relate 9(t) to a M. \ model, with 
e(x) — x, and potential 



a(x) = c J L ,(Xx), 1 < x < 00, 



(166) 



where J„ is the Bessel function of order v, and c and A are 
parameters to be fixed later on. The Mellin transform of 
(I166P (i.e. Fourier for a(q)) yields, 
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o r f u+it ) 

- - MX) = (A/2)-'* ) 2 ( = ±1. 

c v y w ; r(^) 

(172) 

This equation coincides with (|165[) , if we take the minus 
sign in the RHS in (p~72|) and choose 



A = 2tt, c 



(173) 



FIG. 7: Comparison between the absolute values of f(l/2+it) 
and cos B(t). Observe the proximity of the points where both 
quantities vanish. 



Since J\/2{x) — yj2/itx sin (a;), the corresponding poten- 
tial is 



a(x) 



2 sin(27ra;) 



(174) 



a(t) = /°° dx c x~ l+it MXx) 



(167) 



= c 2- 1+lt \- lt 



r(i+^) 

t)r(\+u) ^T~> 1 + 2 ) 1 - 



c_2 



where iF 2 is a hypergeometric function of type (1,2)^2. 
By the Titchmarsh theorem, a(t) is an analytic function 
in the upper half plane. Indeed, the poles of the gamma 
function in the numerator of the first term are cancelled 
out by the poles of the second term. In the limit where 
\t\ — > oo, one gets 



The Jost function T\ it) is found numerically using equa- 
tion eq. (|144p . which involves the Hilbert transform of 
|/i(i)| 2 . Fig. M shows the Argand plot of T\{t), which 
consists of a series of loops passing very close to the 
origin, at those values of t accurately approximated by 
eq. (|165p . For \t\ >> 1 the loops are circles of radius 2 
centered at x — 2. This numerical result can be obtained 
analytically. After a long computation one finds that 



F 1 {t) = l-ce 2i6{ » + 1-+OP), (175) 



which for c = —2 becomes, 



„ (v + it i v + tt , A 2 

hm iF 2 ;H , l + w. 

|i|-oo V 2 ' 2 ' ' 4 



o-Pi 1 + i/, 



A 2 



which is related to the Bessel function 



1 



^(t) =4cos0(t) e lHW +0(- 



(176) 



(168) 

It is interesting to compare (|176[) with (|162[) . which we 
write as 



Mz) 



r(i + i/) 



o-Fi 1 



thus 



r(^) 



and consequently 



r (idd*"\ 
(A/2)- i4 ^^( + J,(A) 



(169) 



|t| » 1, 

(170) 



(xm- u ^^+Mx) 



/l( f)=l + |a(i)~l-£ 

(171) 

A necessary condition for !F\(t) to vanish is that fi(t) 
vanish as well, which in the limit \t\ » 1 is guaranteed 
by 



C(l/2- it) = Z{t) e 



i0(t) 



(177) 



Up to \ jt terms, the phase factor is the same, while Z(t) 
is replaced by cos 9(t), which is precisely the approxima- 
tion that reproduces the smooth part of the zeros (|165[) . 
Hence, to leading order in l/t, the Jost function (|176p can 
be considered as the smooth approximation to £(1/2— it). 



B. Relation to the Berry-Keating regularization 

The Bessel function J v {x) satisfies the second order 
difercntial equation 



which can be rewritten as 
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Re Fi 




x p = E 



FIG. 8: Plot of the numerical evaluation of Ti(t) for the 
potential (|174|) in the range < t < 30. 



FIG. 9: Left: allowed phase space region in the BK semi- 
classical quantization. Right: the dotted line represents the 
classical version of the state f/> a (|182[) . 



(x p) 2 J v (Xx) = h 2 (X 2 x 2 - v 2 )J v (Xx), 
where p = —ih-jn.- From the definition (|33p . 



cq. Q179[) turns into 



cJ v (\x) c sin(27rx) 

y/x 1TX 



(179) 



(180) 



Hi ^ a (x) = K 2 (X 2 x 2 ~v A )^ a (x) 



(181) 



where H — ^fxp^fx is the BKC Hamiltonian 0. Drop- 
ping %j) a in both sides and replacing Hq by xp, one obtains 
a classical version of (|18ip . 



(xp) 2 = ft 2 (A 2 x 2 



±h\< l 



(Ax) 2 



(182) 



which describes a curve in phase space that approaches 
asymptotically the lines p = ±HX. We shall identify these 
asymptotes with the BK boundary in the allowed mo- 
menta \p\ — l p . (see fig. |9|). Recall on the other hand the 
boundary condition x > l x = 1, which combined with the 
previous identification reproduces the Planck cell quan- 
tization condition, 



lp — h\, l x — 1 



l p l x = h\ = 2~nh 



(183) 



where we used that A = 2ir. This interpretation of the 
state ifjaix) shows that the BK boundary |p| — l p is re- 
alized in our model in a dynamical way and not as a 
constraint in phase space. 

To complete this picture, let us try to understand the 
physical meaning of the state tpb for b(x) = 1. Writing 



ip b (x) = sign(x - 1), 
V x 



(184) 



one finds that its time evolution under H is given by 



t) = e~ ltHo Mx) = -7= sign(x - e*), (185) 
which is a kink state associated to the classical trayectory 



x(t) = x e\ p(t) = 0, E = x(t)p(t) = (186) 

Hq does not have a zero energy eigenstate, but based 
on these results one can think of ipi, heuristically as that 
state. The p — line in phase space is the classical ana- 
logue of the state ipb, just like the lines p = ±l p are the 
classical analogue of ijj a . This interpretation allows us to 
understand heuristically, the BK-regularization. Indeed, 
take a particle with energy _E, which at the initial time 
t = is at x = 1 and p = E. Following the classical 
trayectories @ this particle will reach at time t\, the 
boundary p(t\) — l p and then suddenly loose all its mo- 
menta, p(t\ + e) = (see fig. [9]). The phase space area 
involved in this evolution agrees basically with the calcu- 
lation made by Berry and Keating. If this area is a inte- 
ger multiple of 2irh, then there is a bound state. In the 
previous argument, one should strictely reverse the time 
arrow since the classical trayectories generated by 1/ '(xp) 
are the time-reversed of eqs. ((2j) , but the result does not 
change. At the quantum level the existence of bound 
states is due to an interference effect. In the absence 
of this interference the boundary at p = ±l p behaves 
as a "transparent" wall, and the particles do not return 
to their initial position. This situation corresponds to 
Connes picture where all the eigenstates are delocalized. 
In this manner, the Berry-Keating and Connes pictures 
may coexist in a coherent picture both semiclassically 
and quant ally. 

Finally, we would like to make a comment concern- 
ing the wave function (JT5UJ) , and its relation to the von 
Neumann and Wigner potentials, mentioned in section 
^33,38,39,40,41 ^ comrrLOrL feature of these potentials 



is their asymptotic behaviour, sin(r)/r, where r >> 1 
is the radius. The potentials having a positive energy 
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eigenstate form a submanifold. So that a fine-tuning of 
couplings is required. It is interesting to observe the sim- 
ilarity with the potential (|180|) . and the sensitivity to the 
choice of couplings in order to have bound states. 



C. Potentials for the Riemann zeros 

In reference^ Berry and Keating tried to replace the 
semiclassical regularization of xp with quantum bound- 
ary conditions that would generate a discrete spectrum. 
A proposal is to use the dilation symmetry of xp, i.e. 



K x, p^p/K, 



(187) 



where K corresponds to an evolution after time log if 
as indicated by eq.([2]) (see^ for a discussion of the sym- 
metries of H = xp). The Hamiltonian Ho ([6]), is the 
generator of the scale transformations: 
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FIG. 10: Numerical values of Im Tz(E) using two methods: 

1) integrating ()197p in the interval (0, d) with d = 400, and 

2) Summing the series (Tl98)l up to M = 5000. 



ip(Kx) 



1 



ip(x) 



(188) 



Based on this, BK considered a linear superposition of 
the wave function iPe{x) — Cx~ x / 2+lE with the ones 
obtained by integer dilations K = m, 



C 



oo 1 

V - 

jL^i ml/2 



(189) 

C 1 

/2-iE ~ T l/2-iE ^' o _ 



so that the vanishing of (| 1 89[) could be interpreted as 
an eigencondition. However, there is no justification for 
that condition, nor it is clear its physical or geometrical 
meaning. The approach we have been following so far 
is to implement the boundary conditions in a dynamical 
manner, hence it is more natural to impose the symmetry 
under discrete dilations, not on the eigenfunctions but on 
the potentials. Consider the linear superposition of (|180p . 



^2 i>a(mx) 
-l 

oo 

E 



(190) 



m— 1 

oo 



m— 1 



sin(27rmx) c 
rm x 



x] — x + 



where [x] denotes the integer part of x. We have used 
the Fourier decomposition of the sawtooth function [x] — 
x + 1/2. The Mellin transform of the potential a(x) — 
i/xip a (x) associated to (|190[) is given by 



where we have usedi 



C( s ) = s / dx 



\x — X + 



' i Re s > 0. 



s-1 2' 



(192) 

Choosing b(x) = 1, the Jost function Fi(t), associated to 
(I190p . in leading order in t, is given by 



(l+^J - c ( 1 + l) C(l/2-i*)-c 2 5 ao , ao (t), 

(193) 

where S aa>aa denotes the function (|137[) with / = ~g = arj, 
and 



a (t) 



C(l/2- it) 



it 



(194) 



If c = —4, the Jost function T\ is given asymptotically 

by Sao,ao 1 



v ( ~f\ 1 a E2 z ^ 2 p r dt 1 Et z ^ 2 



A further approximation of (|195p is 



(195) 



? 1 (E)/4~Fz{E)= Z{Ef-iP r (:m«i) 

J-co t t-A 



d(t) = / efa;ar i+ " a(x) (191) 



it 



or using that Z(t) = Z(-t) 



T Z {E) = Z{E) 2 -2iE P I — 

Jo n 



dt Z{tf 



t 2 -E 2 ' 



(197) 
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FIG. 11: Numerical values of Im Tz(E) and Im £(1/2 - it) 
in the interval E G (50, 75). 



The Cauchy integral, giving the imaginary part of 
Tz(E), is convergent thanks to the asymptotic behaviour 
\Z(t)\ ~ |t| x /4+e ( e > o), on the critical line ? . The 
integral (|197[) can be computed numerically in the in- 
terval (0, d), with d sufficiently large. The results con- 
verge rather slowly with d. Fig. [TOl shows Im Tz(E) for 
d = 400 in the interval E £ (10, 50). 

An alternative method to find Im Tz(E) is to Hilbcrt 
transform Z{t) 2 , using the well known series expansion 
of the zeta function, 



cw = 21 _ a ~ t Yl 



(-i) n 



Re s > 0. (198) 



71=1 



After a long calculation one finds, 



lmTz(t) = — ilirriM- 



(199) 



1 v> M 
p(t)p(-t) l^n, 

I 

P(*)p(-t) 



M (— l) n + m sign(n— m) 



p(t)-p(-t) ipM 1 



y-»M 2 (- 



/ 2 (l/2 + it){ 



log 2 (n/m)} 



w 



where p(t) — 2 1 / 2+lt — \ and {log 2 (n/m)} is the fractional 
part of log 2 (n/m). FigflOl shows the values of Im Tz{E) 
computed with ea. (|199|) . for M = 5000 in the interval 
E £ (10, 50), which agrees reasonable well with the result 
obtained with the truncated integral (|197[) . For larger 
values of E it is more convenient to use the series expan- 
sion (|199[) . The complete expression of Tz (t) is obtained 
adding to (fT9"9"|) , its real part Z(t) 2 7 



FIG. 12: Argand plane representation of Tz(E) (left) and 
£(1/2 - iE) (right) in the region E SE (50, 75). The loops as- 
sociated to Tz (E) do not generally passed through the origin 
and the are slightly displaced downward. 



T z (t) = limM^o 

J y^M (-] 



(200) 



p7JTpF*T 



,1/2 + it 



E 



M 

n>m— 1 



2 (■ 



I l+p(t)-p(-t) 
+ p(t)p(-t) 
^•.n+rn (^{ 1/2 + it) {log 2 (n/m) } 



Z^n— 1 n 



-(t-^-t) 



As expected, eq. (|200[) does not have poles in the upper 
half plane. The poles arising from the terms proportional 
to l/p(t) cancell each other. 

The real part of TziE) vanishes at the zeros of Z(E). 
The question is wether its imaginary part, given by (|199|l . 
also does. The answer to this question is negative in gen- 
eral, as can be seen from fig. [Til which plots the values of 
Im FziE), and those of 2 Im£(l/2 — iE) in the interval 
(50, 75). Observe that the shape of the two curves is sim- 
ilar, but their zeros do not coincide. Curiously enough, 
their maxima and minima are much closer. Fig. [12] dis- 
plays the Argand plot of Tz{E) and C(l/2 — iE) in the 
interval E 6 (50,75). Observe again the similarity be- 
tween their loop structures. We know of no reason why 
Im Tz (E) should vanish, even asymptotically, at the ze- 
ros of Im C(l/2 — iE) or Z(t). If that were the case, 
then the Riemann zeros would become resonances with 
a life-time increasing asymptotically, but this seems un- 
likely. Our conclusion is that ip a (x) (| 1 90|) is not enough 
to yield the Riemann zeros in the spectra, and that one 
needs a non trivial potential b(x). A concrete proposal 
is to look for a A^2 model yielding a Jost function T{t) 
proportional to C,h (1/2 — it) ( recall (|160|) ). We do not 
see at the moment any obstruction for this realization, 
but one needs additional insights. 

We shall end this section with some comments and 
suggestions. 
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The function Tz{t) (|200p . reminds a two- variable 
version of the zeta function first proposed by Eulcr 



C(si,s 2 ) 



"1>«2>0 



(201) 



if one chooses s\ = 1/2 — it and s 2 = 1/2 + it. 
Eq. (|20ip , together with its multivariable extension, 
called Euler-Zagier zeta functions, have attracted 
much attention in various fields, as knot theory, 
perturbative quantum field theory, etc (see^r^ and 
references therein). The function C{ s i, s 2) satisfies 
the so called shuffle relation 



81,8 2 ) + C(f2, Si) = C(S1)C(S 2 ) " C(*l + 82), (202) 

which amounts in our case to the condition J-z(t) + 
T z {-t) = 2Z{t) 2 . The two variable Euler-Zagier 
function can be meromorphically continued to (C 2 , 
except at the singularities si = 1 and s± + s 2 = 
2, 1, 0, —2, —4, .... Hence, the identification si. 2 = 
1/2 =p it is singular and a proper definition of 
C(si, s 2 ) requires a renormalization, probably along 
the lines of reference^ using Hopf algebras. 

The results obtained in this section can be gen- 
eralized to Dirichlet L-functions with real charac- 
ters ^Sl. We summarize briefly the main results. 
From the functional relation satisfied by the L- 
functions, it follows that the potentials, reproduc- 
ing asymptotically the smooth positions of the ze- 
ros of L(s,x), are given by 



a(x) oc 



sia.(Xx)/y/x, x '■ even 
cos(Ax)/ v / x, x '■ odd 



(203) 



which corresponds to the Bessel functions J v (\x) 
with v = ±1/2. The value A is related to the period 
of the character / as 



A= j, X(n + /)=x(n), Vn e N- (204) 



The generalization of (|190p is 



^— x ' m cos(Amx) 

m=l m— 1 k v y 



(205) 



so that the characters x(m) must be real. The cor- 
responding Mellin transform of a(x) — x 1 / 2 ip a (x) is 
proportional to the Dirichlet function L(l/2—it, x)- 
As in the case of the zeta function, one also needs 
non trivial b potentials to relate the Dirichlet func- 
tions to the Jost functions of the M.2 model. 
VII. CONCLUSIONS 

In this paper we have proposed a possible realization 
of the Hilbert-Polya conjecture in terms of a Hamiltonian 
given by a perturbation of H = xp, or rather its inverse 
X/(xp), by means of an antisymmetric matrix parame- 
terize by two potentials a(x) and b(x). The Schrodinger 
equation can be reduced to a first order differential equa- 
tion, suplemented with boundary conditions, which are 
exactly solvable in terms of a Jost function. In this re- 
spect, our approach is essentially different from the sec- 
ond order approaches to the RH based on standard QM. 

The generic spectrum consists in a continuum of 
eigenenergies which may contain a point like spectra em- 
bbeded in it. We have studied a variety of examples show- 
ing that the existence of a point like spectrum depends 
"critically" on the values of the coupling constants of the 
model. We have found the potentials whose resonances 
approach the smooth Riemann zeros asymptotically. In 
the classical limit these potentials reproduce the Berry- 
Keating semiclassical regularization of H — xp. Imple- 
menting a discrete dilation symmetry on the previous 
potentials, we have obtained a Jost function which, in 
the asymptotic limit, resembles the two- variable Euler- 
Zagier zeta function £(si, s 2 ) with si, 2 = 1/2 =p it. The 
real part of this Jost function vanishes at the Riemann 
zeros but not necessarily its imaginary part. These re- 
sults were derived from a trivial potential b(x) — 1, and 
they suggest that a non trivial choice of b(x), could yield 
a Jost function directly related to the zeta function. A 
natural candidate is ^-((s), with s = 1/2 — it, which 
has the correct analiticity properties. If these potentials 
do exist, then the Riemann zeros would become bound 
states of the model and the RH would follow automat- 
ically. This is how the Hilbert-Polya conjecture would 
come true in our approach. 
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